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Abstract 


A new chronology for the late Miocene of the Mediterranean is presented by combining magnetostratigraphic, 
biostratigraphic (planktonic foraminifera and dinoflagellates) and cyclo-stratigraphic data. Long and continuous upper 
Miocene sections on Gavdos (Metochia section) and Sicily (Gibliscemi section) display cyclic alternations of homogeneous 
marls and sapropels and can be correlated on the basis of their distinct cyclic patterns. The Metochia section yields a good 
paleomagnetic signal and the position of 17 polarity reversals can be determined. The resulting polarity sequence allows an 
unambiguous correlation to the geomagnetic polarity time scale (GPTS). The paleomagnetic signal in the Gibliscemi section 
is too weak to determine a reliable polarity sequence, except for the lowermost part of the section. Detailed biostratigraphic 
analysis results in the identification of 13 planktonic foraminiferal and 9 dinoflagellate bioevents, which can all be accurately 
dated. The Tortonian/ Messinian boundary, defined by the First Regular Occurrence (FRO) of the G. conomiozea group, is 
determined in chron C3Br.ir with an age of 7.12 Ma, according to the GPTS (CK95) of Cande and Kent [1]. 


1. Introduction 


The chronology of the Mediterranean Pliocene 
and Pleistocene has improved substantially during 
the last decade since several integrated magne- 
tostratigraphic and biostratigraphic studies have been 
published on land sections and deep sea cores in the 
eastern Mediterranean [2-6]. All these sections re- 
veal characteristic sedimentary cycle patterns, which 
were shown to be related to the Earth's orbital cycles 
of precession and eccentricity. These sections were 
used to construct an astronomically calibrated polar- 
ity time scale (APTS) for the Plio- and Pleistocene 
[7,8]. Such an APTS proved to be more accurate and 


to have a higher resolution than conventional time 
scales. 

The chronology of the Mediterranean late 
Miocene, however, is still very poor and restricted to 
sediments of late Tortonian-early Messinian age. 
Suitable land-based sections consisting of open ma- 
rine sediments are rare. On Crete, good paleomag- 
netic results were obtained from various sections 
consisting of cyclically bedded marine marls [9,10]. 
These sections, displaying bipartite sedimentary cy- 
cles, do not extend further than 7.7 Ma [11]. Paleo- 
magnetic results from late Miocene sequences in 
Sicily and northern Italy are disappointing. Sections 
represent too short an interval [9], have unsuitable 
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paleomagnetic properties [12], or are completely 
overprinted by a present-day normal field component 
[13]. Finally, magnetostratigraphic studies from ODP 
site 654 in the Tyrrhenian sea gave reasonable re- 
sults but again here the succession is restricted to 
sediments younger than 7 Ma [4]. 

Our present research aims at the development of a 
reliable and accurate chronostratigraphic framework 
for the late Miocene and the extension of the Plio- 
Pleistocene APTS into the Miocene. We have found 
long and continuous upper Miocene sections consist- 
ing of an alternation of marine marls and sapropels 
on the island of Gavdos (Greece) and on Sicily 
(Italy). These sections have been sampled in detail, 
within the framework of the EU-sponsored MIOMAR 
project, and appeared to be very suitable as Mediter- 
ranean reference sections for the late Miocene. 

In this paper, we present a new and detailed 
integrated magnetostratigraphy, biostratigraphy and 
cyclostratigraphy for the upper Miocene from these 
sections on Gavdos and Sicily. Furthermore, we 
incorporate earlier results from sections on Crete for 
the interval straddling the Tortonian/ Messinian 
boundary [11]. The direct calibration of the sedimen- 
tary cycles to the astronomical records and the result- 
ing numerical ages for the individual polarity rever- 
sals and biostratigraphic datum planes is presented in 
Hilgen et al. [14]. 


2. Sections and sampling 


The two main sections used in this study are the 
Metochia section on Gavdos and the Gibliscemi 
section on Sicily which both consist of open marine, 
cyclically bedded marl sequences (Fig. 1). Addi- 
tional data are supplied by sections from Crete 
(Faneromeni, Kastelli) which have previously been 
studied [9—11]. In all sections, sampling was aimed 
at a minimum of four levels per sedimentary (preces- 
sion-related) cycle which corresponds to a resolution 
of approximately 5 kyr. At each level we took four 
cores; two standard-oriented paleomagnetic cores and 
two additional cores for biostratigraphical, sedimen- 
tological and geochemical purposes. At some levels, 
mainly in the laminated beds, it was impossible to 
drill an oriented core and we only took a non-ori- 
ented handsample for biostratigraphic analysis. 


Eea 
jee 
(bees 


AS 


Fig. 1. Location map of the Metochia (MET), Gibliscemi (G/B), 
Kastelli (KA) and Faneromeni (FA) sections (top); detailed maps 
showing the trajectories of the subsections Metochia B and C 
(middle) and Gibliscemi A, B and C (bottom). 


2.1. Gavdos: Metochia section 


The island of Gavdos is located approximately 30 
km south of Crete and represents the southernmost 
emerged part of the Hellenic arc system. The 
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Metochia section is located on the northern part of 
the island and is composed of two subsections 
(Metochia B and C) which are separated by a normal 
fault (Fig. 1; [15]). These subsections are correlated 
on the basis of characteristic lithology patterns. 

The basal part of Metochia B is formed by a 
paleosol which is overlain by sediments character- 
istic of an estuarine environment and shallow-marine 
sands rich in Heterostegina, echinoderma and mol- 
luscs. These so-called Heterostegina sands rapidly 
pass into deep-marine (> 200 m) sediments consist- 
ing of an alternation of hemipelagic marls and 
brownish laminated beds (sapropels) which are num- 
bered M1 to M96 (Fig. 2). In the lower part of 
Metochia B (M7-M11) five successive turbiditic 
sequences occur within or entirely replace the sapro- 
pelitic intervals [16]. In the part of the section below 
the turbidites we succeeded to drill in fresh (blue) 
coloured sediments after removing the weathered 
surface. The middle part of the Metochia B section is 
thoroughly weathered -many samples are charac- 
terised by rusty spots - and even after serious efforts 
no fresh samples could be drilled. 

Metochia C predominantly consists of a rhythmi- 
cally layered succession of hemipelagic marls and 
sapropels. In the upper part of the section, the thick- 
ness of the sedimentary cycles decreases although 
several thin turbiditic beds are intercalated. Metochia 
C is sampled from cycle M36 to M96, resulting in an 
overlap of 20 cycles with Metochia B. The sediment 
is less weathered and brown rusty spots are absent. 
The sapropel-bearing succession in the Metochia sec- 
tion is overlain in stratigraphic continuity by cycli- 
cally bedded diatomites of Messinian age. These 
diatomites, which contain some 38 cycles, are not 
included in the present study. 

We sampled 560 levels over a stratigraphic inter- 
val of approximately 100 m, which corresponds to an 
average spacing of 20 cm. A second sampling trip 
was organised to confirm several important polarity 
reversals and to resample the weathered part of 
Metochia B for magnetostratigraphic purposes only. 


2.2. Sicily: Gibliscemi section 


The Gibliscemi section is located on the southern 
slope of Monte Gibliscemi on the southern part of 


the island of Sicily (Fig. 1). Here, upper Miocene 
sequences start with a cyclic alternation of homoge- 
neous hemipelagic marls and sapropels, followed by 
the diatomites of the Tripoli Formation (Messinian) 
and capped by the Calcare di Base and evaporites 
(gypsum) of the Gessoso-Solfifera Formation. 

The basal part of the Gibliscemi section is tectoni- 
cally disturbed and overlies the so-called ‘‘argille— 
scagliose’’ which forms the decollement level of 
thrust nappes. Various shear planes cause important 
deformations/hiatuses and no continuous section 
could be sampled along a single trajectory. We 
succeeded in obtaining a relatively undisturbed and 
continuous succession -in which the sapropels are 
numbered (G1—G99, Fig. 3) -by sampling three 
separate subsections (Gibliscemi A, B and C; Fig. 1) 
which can be correlated on the basis of their distinct 
characteristic cyclic patterns. In our composite sec- 
tion, deformation is only bedding-parallel and (pre- 
dominantly) restricted to sapropels. Only in the top 
part of the section, a shear plane is found which runs 
obliquely to the bedding, in the interval preceding 
the diatomites of the overlying Tripoli Formation. 

A. total of 620 levels has been sampled in the 
three subsections, corresponding to an average spac- 
ing of 20 cm. A second sampling trip was made to 
focus on the lowest part of the section. We resam- 
pled the interval comprising cycles G1—-G25 in the 
Gibliscemi A and C sections to obtain a more reli- 
able magnetostratigraphy. 


3. Cyclostratigraphy 


All sections consist of open-marine sediments that 
show cyclic alternations of either whitish-coloured 
carbonate-rich and grey-coloured carbonate-poor 
marls (Faneromeni section; lower part), or of homo- 
geneous marls and brownish-coloured laminated beds 
termed sapropels (see [11]). Especially these sapro- 
pels reveal characteristic patterns - with both small- 
scale and large-scale clusters in addition to the indi- 
vidual sapropels - which allows the sections to be 
correlated (cyclostratigraphically) in detail. The sedi- 
mentary cycle patterns of the studied sections and 
their correlation to the astronomical target curves are 
discussed in detail by Hilgen et al. [14]. 
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4. Magnetostratigraphy 


Thermal demagnetisations were performed in a 
magnetically shielded, laboratory built furnace. The 
natural remanent magnetisation (NRM) was mea- 
sured on a 2G Enterprises cryogenic magnetometer. 
At least one specimen per sampling level was ther- 
mally demagnetised using temperature increments of 
30—50°C. Furthermore, we performed some rock 
magnetic tests to identify the dominant carriers of 
the magnetism, including acquisition of an isother- 
mal remanent magnetisation (IRM) and subsequent 
demagnetisation of this IRM. The IRM was induced 
in a pulse magnetiser and measured on a digitised 
spinner magnetometer based on a Jelinek JR3 driver 
unit. After each step we measured the low field 
susceptibility which was measured on a Kappabridge 
KLY-2. 


4.1. Gavdos: Metochia section 


Thermal demagnetisation of the samples from the 
Metochia section revealed different types of demag- 
netisation behaviour, which appeared to be rather 
constant throughout specific intervals. Viscous be- 
haviour, when heated above temperatures of 400°c, 
is shown by the marly interval below the turbidites 
(M1—M6). Good-quality demagnetisation (linear de- 
cay) is observed in the Heterostegina sands, the 
turbidites and four following cycles (M7—M15), the 
upper part of Metochia B (M43—M56) and the lower 
part of Metochia C (M36-M56, with exception of 
M49). Low-intensity clusters at temperatures be- 
tween 240 and 600°c are observed in the middle part 


of Metochia B (M16-M42) and the upper part of 
Metochia C (M57-M93) but show different be- 
haviour in these two subsections. We performed 
some rockmagnetic tests on a number of samples 
from each interval to characterise the magnetic prop- 
erties and to distinguish primary from secondary 
components. On the basis of these experiments, we 
divided the samples in four groups (Fig. 4). For each 
group, all selected samples show the same typical 
behaviour. Therefore, we assume that all samples 
showing the same demagnetisation behaviour can be 
attributed to the same group. 


4.1.1. Group A 

Group A is characterised by the dominance of a 
low-coercivity mineral with maximum blocking tem- 
peratures below 600?C (Fig. 4a) indicating the pres- 
ence of magnetite. Heating of the samples to temper- 
atures above 390°C results in a sudden increase of 
susceptibility which continues up to temperatures of 
500°C. This behaviour was observed earlier in sub- 
oxic to anoxic sediments and may be the result of 
oxidation of an iron sulphide like pyrite [17]. The 
decrease in susceptibility between 500 and 600°C 
can be ascribed to oxidation of magnetite to hematite. 

Samples attributed to group A are found in the 
lowermost 6 cycles (M1—M6) of Metochia B, be- 
tween the Heterostegina sands and the turbidites. 
NRM intensities range from 0.04 to 2 mA /m, initial 
susceptibilities from 90 to 180 x 1076 SI. The char- 
acteristic remanent magnetisation (ChRM) is usually 
largely removed at temperatures ranging from 300 to 
390°C. Demagnetisation diagrams for this interval 
show both normal and reversed polarities (Fig. 5c,d) 


Fig. 2. Polarity zones, lithology, cycle numbers and ranges of planktonic foraminifera, Bolboforma and dinoflagellates of the Metochia 
section. In the polarity column black (white) denotes normal (reversed) polarity interval and shaded interval denotes zone of undefined (only 
secondary components) polarity. Closed (open) symbols denote reliable (unreliable) directions. v = group A samples; @ = group B 
samples; 9 = group C samples; A = group D samples; * = directions obtained by applying the great-circle method [18]. Lithology column 
shows cyclic alternations of homogeneous marls (white) and sapropels (black); shaded interval represents turbiditic intercalations [16]. Right 
part of the figure shows biostratigraphic data given in terms of absence /presence. Numbered bioevents (encircled) refer to (1) highest 
regular occurrence (hro) of dextral N. acostaensis, (2) lowest common occurrence (Ico) of G. menardii 4, (3) Last Common Occurence 
(LCO) of large-sized C. parvulus, (4) frequency shift (fs) in C. parvulus, (5) lowest regular occurrence (lro) of S. seminulina, (6) hro of S. 
seminulina, (7) LCO of G. menardii 4, (8) LO of C. parvulus, (9) First Occurrence of G. menardii 5 (FO), (10) First Regular Occurrence 
(FRO) of the G. conomiozea group, and (11) FO of G. nicolae, (I) FO of the G. conomiozea group. Arrows indicate positions of bioevents. 


+ = large-sized types of C. parvulus. 
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suggesting a primary origin for this component. Fur- 
ther demagnetisation at higher temperatures results 
in a mainly randomly directed viscous component. 


4.1.2. Group B 

IRM acquisition of group B samples indicate the 
presence of both a low-coercivity (magnetite) and a 
high-coercivity mineral (Fig. 4b). IRM demagnetisa- 
tion shows that at least one component persisted at 
temperatures higher than 600°C. Approximately at 
680°C the IRM is totally removed which indicates 
that the high-coercivity mineral is hematite. Most 
samples show no significant change in susceptibility 
during the heating experiments. Samples which are 
attributed to group B are observed in the Heteroste- 
gina sands, cycles M7—15 and M43-56 in Metochia 
B and in cycles M36-MS6 (except M49) in Metochia 
C. 

NRM intensities range from 1 to 7 mA/m, sus- 
ceptibilities from 100 to 200 x 107% SI. Thermal 
demagnetisation diagrams show that in most samples 
a large secondary present-day field component is 
removed at 240°C (Fig. 5a,j,k). Demagnetisation at 
higher temperatures reveals both normal and re- 
versed components and again suggests a primary 
origin of the ChRM. Demagnetisation diagrams are 
generally of good quality and in most cases the 
polarity can reliably be determined. Some samples 
do not show a linear decay to the origin but move 
successively toward the reversed quadrant of the 
Zijderveld diagram (Fig. 5e). We assume that the 
primary component is partly overprinted by a sec- 
ondary component caused by overlapping blocking 
temperature spectra. When plotted on an equal-area 
diagram, remanence vectors progressed along a great 
circle toward a southerly upward direction indicating 
removal of a normal phase from a reversed primary 


component. In these cases a best-fitting great-circle 
plane was used to estimate the characteristic direc- 
tions according to the method of McFadden and 
McElhinny [18] (Fig. 5f. 


4.1.3. Group C 

IRM acquisition and demagnetisation of group C 
samples show that the magnetisation is here domi- 
nated by a high-coercivity mineral (Fig. 4c). IRM- 
acquisition curves clearly differ from the group A 
and B curves. The blocking temperatures range from 
660 to 680°C which is characteristic for hematite. No 
major changes in susceptibility are observed during 
heating up to 700°C. Samples of group C are found 
in cycles M16-M43 in Metochia B. 

NRM intensities range from 1 to 3 mA /m, initial 
susceptibilities from 80 to 160 x 10 5 SI. Demag- 
netisation diagrams show that a large normal compo- 
nent is removed at temperatures of 240°C. At higher 
temperatures a low-intensity cluster is observed, al- 
ways indicating normal magnetisation directions (Fig. 
5h,i). The samples from cycles M40 and M43 in 
Metochia C, however, give clearly reversed primary 
directions whereas demagnetisation of the samples 
from the same cycles in Metochia B results in clus- 
ters which are characteristic for group C (Fig. 5i,j). 
This indicates that the high-coercivity mineral 
(hematite) is of secondary origin and most likely 
generated by subrecent weathering of the sediment. 


4.1.4. Group D 

Samples attributed to group D again show a low 
coercivity mineral as the dominant carrier of the 
remanence (Fig. 4d). Maximum blocking tempera- 
tures of approximately 580°C point to the presence 
of magnetite. The difference with group A is the 
response of susceptibility during heating of the sedi- 


Fig. 3. Polarity zones, lithology, cycle numbers and ranges of planktonic foraminifera, Bolboforma and dinoflagellates of the Gibliscemi 
section. Biostratigraphic data of planktonic foraminiferal marker species are based on surveying a standard number of fields (27 out of 45) 
on a rectangular picking tray and semi-quantitatively presented in terms of absence, trace (< 3 specimens per 9 fields of picking tray), rare 
(3-10), common (10-30) and frequent (> 30) indicated by increasing bar thickness. For numbered foraminiferal bioevents see caption to 
Fig. 2. Semi-quantitative distribution of dinoflagellate marker species is based on 200-counts and presented in terms of absence, trace 
(1-5%), rare (5-10), common (10-25) and frequent (> 25) indicated by increasing bar thickness. Encircled dinoflagellate events refer to 
(a) LO of C. sp. A, (b) LCO of M. sp. A, (c) FO of I. maculatum, (d) LO of P. golzowenze, (e) FO of C. labradorii, (f) LO of M. sp. A, 
(g) of LO C. labradorii, (h) LCO of T. psilatum, and (i) LO of L. truncatum. 
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Fig. 4. Examples of rock magnetic tests for selected samples of the four groups (A, B, C and D) of the Metochia section. In the IRM 
acquisition diagrams, the initial steep rise (0-300 mT) points to magnetite, the gradual increase at high fields suggests the presence of 
hematite (see especially group C). Stepwise thermal demagnetisation of the normalised NRM also show the presence of magnetite (A and D) 
and hematite (B and C). Low field susceptibilities display chemical reactions during thermal demagnetisation. Note the sudden increase at 
approximately 400°C in group A. 
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ment. In group D, no increase in low field suscepti- 
bilities is observed between temperatures of 390- 
500°C. Samples attributed to group D are found in 
cycle M57 of Metochia B and in cycles M57-M91 
of Metochia C. 

NRM intensities range from 0.17 to 2.5 mA/m, 
susceptibilities from 35 to 150 x 1076 SI. The lower 
values of the initial susceptibility are probably re- 
lated to a higher CaCO, content. Demagnetisation 
diagrams show that a secondary low-temperature 
(LT) component is removed at temperatures of 240°C 
(Fig. 5k—o). In the lower part of the section this LT 
component has a normal present-day field direction, 
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whereas in the upper part of the section this LT 
component has a reversed direction (Fig. 5m-o). 
After removing the LT component, a relatively small 
but well-determined high-temperature (HT) compo- 
nent reveals a consistent pattern of both normal and 
reversed polarities (Fig. 5). We suggest that the HT 
component (> 240°C) represents the primary ChRM 
component and the LT component a secondary 
post-depositional overprint. The reversed polarity of 
the LT component might be generated by late-di- 
agenetic processes, for instance when the sediment 
became uplifted/exposed during a period of re- 
versed polarity. The rock magnetic results indicate 
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Fig. 6. (A) Examples of IRM acquisition, subsequent IRM demagnetisation and low-field susceptibilities for selected samples of the 
Gibliscemi section. (B-D Thermal demagnetisation diagrams. See also captions to Figs. 4 and 5. 
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that the secondary component in Metochia C is a low 
coercivity (probably magnetite) component. 

The ChRM directions and polarity zones of the 
composite Metochia section show that at least seven- 
teen polarity reversals are recorded (Fig. 2; see also 
Appendix A). The magnetostratigraphic results from 
the partly overlapping Metochia B and C sections 
confirm the cyclostratigraphic correlation. No pri- 
mary components could be determined in the middle 
part of the Metochia B section (M16—M42). In this 
part of the section the secondary component totally 
overprints the primary component and hence no orig- 
inal polarity can be determined. The lower part of 
Metochia C is less weathered and yields good results 
for cycles M40—M42. The paleomagnetic signal of 
cycles M36-M39 is too weak to determine a reliable 
polarity. 


4.2. Sicily: Gibliscemi section 


The NRM intensities of the Gibliscemi section are 
very low, 50-150 uA /m and initial susceptibilities 
range from 100 to 200 x 10 — 6 SI. Only for the 
lower part of the section reliable magnetostrati- 
graphic data could be derived (Fig. 3). The best 
demagnetisation diagrams were given by samples of 
the lowermost part (cycles G1—G25) of the section. 

IRM acquisition and subsequent demagnetisation 
for selected samples of the lower part of Gibliscemi 
A indicate the presence of a low-coercivity mineral 
with a maximum unblocking temperature of 560- 
600°C (Fig. 6a). This suggests that magnetite is the 
dominant carrier of the magnetism. Heating of the 
samples results in an increase of the low field sus- 
ceptibility between temperatures of 400 and 480*C 
after which a decrease to 600°C is observed. These 
results generally show the same characteristics as the 
samples of Group A of the Metochia section. 

Thermal demagnetisation often reveals that a sec- 
ondary, normally directed component is removed at a 
temperature of 240°C (Fig. 6c) although this compo- 
nent is not always present (Fig. 6e). The ChRM 
component is subsequently demagnetised at tempera- 
tures of approximately 400°C. Demagnetisation at 
higher temperatures results in a large and randomly 
directed viscous component, compatible with an in- 
crease in susceptibility probably caused by the oxida- 
tion of pyrite to magnetite [17]. Demagnetisation 


diagrams are often difficult to interpret because of 
the very weak NRM (Fig. 6b,f,h). We interpreted a 
direction as reliable only if a linear decay to the 
origin was observed. The reliable ChRM directions 
show an average clockwise rotation; 20^ for the 
normal polarities and 25° for the reversed polarities. 
Clockwise rotations on Sicily have been found ear- 
lier by Scheepers and Langereis [19] who estimated 
an average clockwise rotation of 34° for early 
Pliocene sections from the Caltanissetta basin. 

Four polarity reversals are recorded in the lowest 
part of the section (Fig. 3; see also Appendix A). 
The results from the Gibliscemi A and C sections are 
in agreement with each other and confirm the posi- 
tions of the polarity reversals. Unfortunately, the 
upper two reversals cannot be pinpointed accurately 
in one sedimentary cycle because of the very weak 
paleomagnetic signal. 


5. Biostratigraphy 


Planktonic foraminiferal biostratigraphy played an 
essential role in confirming the cyclostratigraphic 
correlation between Gibliscemi and Metochia since 
magnetostratigraphic and dinoflagellate control is 
poor in one or the other section. Planktonic 
foraminiferal biostratigraphic correlations are based 
on the stratigraphic distribution of six species plus 
the coiling ratio of Neogloboquadrina acostaensis in 


,337 samples from Gibliscemi and 553 from Metochia 


(see Figs. 2 and 3). The distribution of these six 
species is based on surveying a standard number of 
fields (27 out of 45) on a rectangular picking tray 
and specified in semi-quantitative terms (frequent, 
common, rare, trace) for Gibliscemi and in qualita- 
tive terms (presence, absence) for Metochia. Most 
species display an intermittent distribution pattern. 
This holds in particular for the sinistrally-coiled 
Globorotalia menardii 4 in the lower part of both 
sections (G1-G50 and M1-M45). This species oc- 
curs more regularly between G51-G80 and M46- 
M68 but then vanished until its brief re-occurrence 
in G85 of Gibliscemi and in M69 of Metochia. The 
top of this brief re-occurrence event is equated with 
the LCO of G. menardii 4 even though this species 
shows a final and extremely brief influx in the upper 
part of the range of G. menardii 5 [20]. Fig. 2 shows 
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that the sinistral G. menardii 4 is followed in succes- 
sion by the dextral G. menardii 5 and the sinistral 
Globorotalia conomiozea group. This succession (al- 
though differently labeled) has been noticed in many 
locations in the Mediterranean [21—25], and the ab- 
sence of G. menardii 5 in Gibliscemi, therefore, 
points to a major hiatus at 108 m (Fig. 3). The 
unkeeled globorotaliid group lacks systematic 
changes that can be used for correlating both sec- 
tions. An exception is the FO of Globorotalia nico- 
lae -a dextrally-coiled and biconvex form with in- 
flated chambers - which is a useful bioevent for 
correlating lower Messinian sequences in the 
Mediterranean region. The FO of G. nicolae in 
Metochia shortly preceeds the first occurrence of 
diatomites and follows a relatively long interval in 
which unkeeled globorotaliids are absent. The top of 
Gibliscemi reaches into this absence interval. An- 
other useful species is Catapsydrax parvulus, which 
is the last representative of the Catapsydrax group. 
Catapsydrax parvulus is a senior synomym of 
Globorotalia falconarae, which, in the Mediter- 
ranean ranges up into Zone N16 [26,27]. Catapsy- 
drax parvulus is common to frequent in cycles G1- 
G37 of Gibliscemi and in M1-M29 of Metochia and 
its maximum size is significantly larger than in any 
younger level (with the exception of a brief reoccur- 
rence in G72 and M61; see Figs. 2 and 3). The LCO 
of these large-sized C. parvulus types is a useful 
bioevent for the (eastern) Mediterranean. Catapsy- 
drax parvulus is extremely scarce from the LCO 
level of the large-sized types up to cycle G37 and 
M29, but becomes again more frequent from cycles 
G53 and M47 up to its LO level (M71) in Metochia 
(Figs. 2 and 3). 

Sphaeroidinellopsis seminulina is common to fre- 
quent in a discrete interval which ranges from G72 
to G80 and from M61 to M68 (Figs. 2 and 3). The 
base of this interval coincides with the brief re-oc- 
currence of large-sized C. parvulus types, while the 
top corresponds with the base of the longer-term 
absence interval of G. menardii 4 preceding its LCO 


level. Scattered occurrences of S. seminulina are 
found beyond this interval in cycle M10 and M68- 
M69. The last species with biostratigraphic signifi- 
cance is N. acostaensis. Coiling of this species is 
dominantly dextral from G1 to G10 and from MI to 
M6, whereafter coiling changes in a series of high- 
amplitude fluctuations up to G19 in Gibliscemi and 
M15 in Metochia. Above that level dominant sinis- 
tral coiling prevails, although at some levels coiling 
reverts to dominant dextral. 

The stratigraphic distribution of the species listed 
above provides 11 first-order bioevents (summarized 
in Appendix B), which are perfectly isochronous 
compared to the cyclostratigraphic framework and 
which substantially improve the biostratigraphic res- 
olution in this time span. However, many more 
isochronous surfaces can be drawn between both 
sections on the basis of second-order bioevents. Ex- 
amples of such second-order bioevents are several 
short-term incursions of Globoquadrina altispira, 
the many short-term absence intervals of G. menardii 
4, and the brief reoccurrence of large-sized C. parvu- 
lus types. One of the most remarkable second-order 
bioevents is the single short incursion of abundant 
Globoquadrina dehiscens in cycle G32 of Giblis- 
cemi and in the equivalent cycle M28 of Metochia. 

Along with the planktonic foraminiferal analysis 
we recorded the absence and presence of 5 species of 
the calcareous microfossil group of Bolboforma 
(Protophyta, [28]). Figs. 2 and 3 show that the 
sections Metochia and Gibliscemi have three species 
in common but their stratigraphic distribution shows 
little similarity. This is surprising because the oppo- 
site is true for the distribution of the planktonic 
foraminifers. We conclude, therefore, that this group 
has little biostratigraphic significance in this part of 
the Mediterranean and for the time slice studied. 

Dinoflagellate cysts are abundant in the Giblis- 
cemi section, but rare in Metochia. Actually, 75% of 
the samples in Metochia was barren possibly due to 
intense weathering. 200-counts on 72 samples from 
below the hiatus in Gibliscemi yielded 7 species with 


Fig. 7. Stratigraphic framework showing the magnetostratigraphy, integrated biostratigraphy and cyclostratigraphy of the Kastelli, 
Faneromeni, Metochia and Gibliscemi sections and the correlation of the magnetostratigraphy to the geomagnetic polarity time scale of 
Cande and Kent [1]. Cyclostratigraphic correlations are based on sedimentary cycle patterns. 
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(potential) biostratigraphic significance (Fig. 3). The 
LO of Cerebrocysta sp. A in cycle G6 has also been 
recorded by Powell [29] in the lower part of the 
Tortonian type section along the Rio Mazzapiedi— 
Castellania and might be a useful bioevent for (at 
least) the Mediterranean region. The LO of Palaeo- 
cystodinium golzowense in G41 seems to have 
(semi-)global significance since it has been recorded 
in a number of upper Miocene sequences from out- 
side the Mediterranean as well (e.g. [30—33]). Men- 
dicodinium sp. A has its LO in G61. The highest 
occurrence of Labyrinthodinium truncatum in G78 
represent its true last occurrence since this species is 
absent in the Faneromeni section, which covers the 
interval missing in Gibliscemi due to a hiatus. The 
LO of L. truncatum is a useful bioevent since it has 
been recorded in a number of late Miocene sections 
from outside the Mediterranean as well [30—33]. 
Tectatodinium psilatum ranges up into the 
Faneromeni section but its frequency drop in G68 is 
a very distinct and probably useful bioevent. Whether 
the FO of Impagidinium maculatum in G21 and the 
limited range of Corrudinium labradorii (from G58 
to G66) have biostratigraphic significance will be- 
come evident after studying additional late Miocene 
sections from within and outside the Mediterranean. 


6. Integrated stratigraphic framework 


The data presented in this paper provide a high- 
resolution integrated stratigraphic framework for the 
Mediterranean upper Miocene. Since the Gibliscemi 
section recorded a large hiatus in its upper part and 
also did not provide a good magnetostratigraphy, we 
incorporated the earlier data from the Faneromeni 
and Kastelli sections of Crete which showed reliable 
results for the corresponding interval [10,11]. In 
general, the  magnetostratigraphic and  biostrati- 
graphic data are consistent with, and thus essentially 
confirm, the cyclostratigraphic correlations, although 
minor discrepancies exist (Fig. 7). The polarity re- 
versals of the upper part of the Metochia section are 
consistently found at a slightly lower (cyclo)strati- 
graphic position than at Faneromeni and Kastelli, as 
confirmed by the biostratigraphic data. Furthermore, 
the consistent offset of generally 1/2 to 1 cycle 
suggest that this discrepancy is probably not caused 


by a misinterpretation of the paleomagnetic data. 
Diagenetic processes, related to the paleoredox con- 
ditions, may cause a delay of the acquisition of the 
NRM, which is dependent on the lithology of the 
sediment [34]. Such a delayed acquisition is a likely 
explanation for the observed discrepancies between 
Metochia and Crete. 


7. Correlation to the GPTS 


Since the results of all sections are generally in 
good agreement, we can use the polarity sequences 
to construct a chronostratigraphic framework and to 
assign numerical ages to the biostratigraphic events. 
Correlation of the established polarity patterns (Ap- 
pendix A) to the GPTS has been made by reference 
to the most recent published version CK95 [1]. The 
results of the Cretan sections were correlated earlier 
to the older CK92 [35] time scale [11]. In CK95, the 
polarity patterns are identical to CK92; differences 
are caused by the incorporation of new and addi- 
tional calibration points which results in slightly 
older ages for the polarity reversals. Hence, the 
correlation of the polarity patterns of the Faneromeni 
and Kastelli sections to CK95 remains the same but 
only results in a minor change in age (Fig. 7). This 
was also suggested by the modified CK92 MOD 
time scale of Krijgsman et al. [11]. 

In contrast to the Cretan sections, the Metochia 
section could be extended further downward which 
resulted in the recording of eight additional polarity 
reversals (Fig. 7). Correlation to CK95 shows that 
the long normal interval (M57—M68) corresponds to 
chron C4n.2n and the first short normal interval 
(M51-M323) to chron C4r.1n. The second short nor- 
mal interval (M41—M42) cannot be correlated to any 
chron in CK95. It might correspond to the so-called 
cryptochron C4r.2r-1 which has a duration of 16 kyr 
according to CK95. This short polarity interval is 
also found in ODP Leg 138 [36]. The duration of this 
interval, estimated from thickness and sedimentation 
rate, ranges from 33 to 58 kyr in Schneider [36], 
which is in agreement with our results that indicate a 
duration of 30—40 kyr (1.5-2 sedimentary cycles). 
Hence, a duration of 16 kyr in CK95, based on ocean 
floor anomaly patterns, seems too short. Schneider 
[36] also recorded another cryptochron, which he 
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called C4r.1r-1, that is not recorded in the Cande and 
Kent [1,35] time scale. According to our results, only 
reversed polarities are recorded in the corresponding 
time interval. There is one level (M49), however, for 
which we obtained normal polarities, but rock mag- 
netic experiments suggest that we are dealing here 
with a normal (hematite; group B) overprint. This 
level thus reveals the same rock magnetic character- 
istics as the long remagnetised interval (M16—M39) 
in Metochia B. This remagnetised interval is marked 
by the (regular) occurrence of numerous sapropels 
whereas the remagnetised M49 level occurs in a 
small group (M48-M50) of sapropels. In this re- 
spect, it is suspect that two short normal polarity 
intervals coincide with small sapropel groups (i.e. 
M41-M44 and M51-M54). Their rock magnetic 
characteristics, however, suggest a primary origin of 
these normal polarities. 

No reliable data are derived from the middle part 
of the Metochia B section (M16-M39) because the 
primary component is probably totally overprinted 
by a normally directed secondary weathering compo- 
nent. Also in the Gibliscemi section no polarities 
could be determined for the time-equivalent interval. 
The combined results of Gibliscemi and Metochia 
show that four polarity reversals are recorded at the 
base of the sections, which determine two short 
normal polarity intervals. The only likely correlation 
to CK95 of the 20 kyr (M1; G5) and 50 kyr (M14- 
M15; G18-G19) normal polarity intervals is to the 
subchrons C4Ar.1n and C4Ar.2n (Fig. 7). According 
to this correlation, it appears that we did not record 
the upper and lower reversal boundaries of chron 
C4An. Our correlation is furthermore confirmed by 
the duration of the interval between the reversal 
boundaries C4r.1n (o) and C4Ar.1n (y), which is 973 
kyr in CK95 and 1040 kyr based on the estimated 
number of (4-50 precession related) sedimentary 
cycles. 

The established magnetostratigraphic correlation 
results in accurate ages, according to the CK95 time 
scale, for the following bioevents (Appendix B); (1) 
highest regular occurrence (hro) of dexral N. 
acostaensis (chron C4Ar.2r; 9.53 + 0.01 Ma); (2) 
lowest common occurrence (Ico) of G. menardii 4 
coincident with the last of a series of high-amplitude 
changes in 96 dextral N. acostaensis (C4Ar.in (y); 
9.24 + 0.01 Ma); (3) LCO of large sized C. parvulus 


types (8.81 + 0.01 Ma); (4) frequency shift in C. 
parvulus from nearly absent to commonness (C4r.2r; 
8.40 + 0.02 Ma); (5) the Iro of S. seminulina, coinci- 
dent with a brief re-occurrence of large-sized C. 
parvulus types (C4n.2n; 7.85 + 0.02 Ma); (6) the 
highest regular occurrence (hro) of S. seminulina, 
coincident with the temporary disappearance of G. 
menardii 4 preceding its LCO (C4n.2n(y); 7.65 + 
0.03 Ma); (7) the LCO of G. menardii form 4 
(C3Br.3r; 7.40 + 0.02 Ma); (8) the last occurrence 
(LO) of C. parvulus (C3Br.2r; 7.32 + 0.01 Ma); (9) 
the FO of G. menardii form 5 (C3Br.2r; 7.23 + 0.01 
Ma); (10) the First Regular Occurrence (FRO) of the 
G. conomiozea group (C3Br.1r; 7.12 +0.01 Ma), 
(11) the FO of G. nicolae (C3Ar; 6.68 + 0.01 Ma). 
The dinoflagellate bioevents are only determined in 
the Gibliscemi section which could not be directly 
dated since this section did not provide a good 
magnetostratigraphy. Numerical ages for these events 
are calculated after a cyclostratigraphic correlation of 
the polarity reversals in the Metochia section to their 
corresponding position in Gibliscemi (see Appendix 
B). 


8. The age of the T/M boundary 


The level at 6 m below the base of the Tripoli 
Formation in the Falconara section (Sicily) has been 
proposed in 1979 by Colalongo et al. as the Global 
Stratotype Section and Point (GSSP) for the Torto- 
nian/ Messinian (T/M) boundary and serves as such 
already for 16 years, but the boundary has never 
been formally defined. The position of the proposed 
GSSP for the T/M boundary is marked by the FO of 
G. conomiozea [37]. 

Globorotalia conomiozea belongs to a late 
Miocene group of keeled globorotaliids characterized 
by sinistral coiling and a reniform chamber outline in 
spiral view [20]. Relatively flat morphotypes are 
named G. dalii and G. miotumida (with G. dalii 
being a junior synomym of G. miotumida), whereas 
conical forms are labeled G. conoidea, G. 
conomiozea, G. saphoa, and G. mediterranea de- 
pending on the degree of test thickness, axial angu- 
larity, and the number of chambers. All these mor- 
photypes, however, completely intergrade, which -in 
addition to the similarity in chamber outline and 


490 W. Krijgsman et al. / Earth and Planetary Science Letters 136 (1995) 475—494 


coiling -is the reason why they all have been lumped 
into one single category called the G. conomiozea 
group [20] or the G. miotumida group [24]. 

The late Miocene G. conomiozea group lived at 
the mid-to high latitudes of both hemispheres [38] 
and invaded the Mediterranean near the T/M bound- 
ary in response to increased climatic cooling [39,24]. 
This invasion is preceded, however, by an earlier 
short-term invasion of thick-shelled and conical rep- 
resentatives of this group at two closely spaced 
levels slightly above the lowest regular occurrence of 
S. seminulina (Figs. 2 and 3). The G. conomiozea 
group coexisted with G. menardii 4 during the first 
invasion, but was the sole keeled globorotaliid in the 
Mediterranean during the second invasion near the 
T/M boundary. The sudden spreading of the G. 
conomiozea group in the Mediterranean during this 
second invasion (= FRO of the G. conomiozea 
group) is taken by Langereis et al. [10], Sierro et al. 
[24], and Krijgsman et al. [11] to mark the T/M 
boundary. The earliest representatives of the G. 
conomiozea group during this second invasion con- 
sist of the flat G. miotumida type [23-25]. The 
conical types arrived later and that is the reason why 
some authors (e.g. [23,25]) place the T/M boundary 
at the first occurrence level of these conical (— G. 
conomiozea) types in accordance with the original 
proposal of Colalongo et al. [37]. These conical types 
first occur in section Metochia in cycle M73, shortly 
after the FRO of the G. conomiozea group, but their 
regular occurrence starts in cycle M76 in the younger 
part of chron C3Bn. Hodell et al. [25] might have 
equated the latter level with the T/M boundary, 
which would explain why they placed this boundary 
in chron C3Bn, whereas Langereis et al. [10] and 
Krijgsman et al. [11] placed the T/M boundary 
midway in chron C3Br.1r. The earliest representa- 
tives of the G. conomiozea group at the proposed 
T/M boundary level in the Falconara section contain 
already some conical types, but their regular occur- 
rence starts at about 8 m above the proposed bound- 
ary level. This suggests that the proposed T/M 
boundary level in Falconara is associated with a 
small hiatus and, therefore, might not be a suitable 
candidate for defining the T/M boundary. In antici- 
pation of a formal definition we place the T/M 
boundary at the FRO of the G. conomiozea group (in 


accordance with [10] and [11]) because this bioevent 
reflects the most pronounced change in late Miocene 
globorotaliids and is - irrespective of species concept 
and labeling procedures -easy to identify within and 
immediately outside the Mediterranean. The T/M 
boundary defined in this way has an age of 7.12 Ma 
based on linear interpolation between the younger 
end of chron C3Br.1n and the older end of chron 
C3Bn in CK95. 


9. Conclusions 


A high-resolution stratigraphic framework for the 
late Miocene in the (eastern) Mediterranean is estab- 
lished by combining magnetostratigraphic, biostrati- 
graphic and cyclostratigraphic data from marine sec- 
tions on Gavdos, Sicily and Crete. Seventeen polar- 
ity reversals are recorded which allows an unam- 
biguous correlation to CK95 (Fig. 7). This correla- 
tion is used to calculate ages for eleven major plank- 
tonic foraminiferal and nine dinoflagellate bioevents 
(Appendix B). The Tortonian/Messinian boundary, 
placed at the FRO of the G. conomiozea group, is 
determined in chron C3Br.1r and has an age of 7.12 
Ma according to CK95. 
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Stratigraphic levels of polarity reversals recorded in the studied sections. 
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C3An.2n (0) 
C3Bn (y) 
C3Bn (o) 
C3Br.1n (y) 
C3Br.1n (0) 
C3Br.2n (y) 
C3Br.2n (0) 
C4n.1n (y) 
C4n.in (0) 
C4n.2n (y) 
C4n.2n (0) 
C4r.1n (y) 
Cár.1n (0) 
C4r.2r-1 (y) 
C4r.2r-1 (0) 
C4An (y) 
C4An (0) 
C4Ar. ln (y) 
C4Ar.in (0) 
C4Ar.2n (y) 
C4Ar.2n (o) 


57,43 - 57.50 
31.60 - 32.41 
28.63 - 29.73 
25.84 - 26.18 
15.35 - 15.60 
11.80 - 13.40 
8.35 - 9.00 


55.72 - 57.00 
46.83 - 47.25 
43.40 - 43.90 
40.57 - 41.03 
29.19 - 29.67 
26.79 - 27.39 
24.25 - 25.40 
16.75 - 17.39 
13.29 - 13.97 


92.57 - 9321 
90.60 - 90.74 
89.69 - 89.91 
88.72 - 89.03 
8491 - 85.02 
83.89 - 84.00 
83.09 - 83.28 
80.23 - 80.87 
78.83 - 79.87 
68.86 - 69.26 
62.99 - 63.90 
61.36 - 61.82 
50.79 - 51.60 
49.29 - 49.90 


27.35 - 27.41 
5.63 - 573 
2.55 - 3.10 


17.99 - 20.46 
12.09 - 14.25 
4.55 - 5.05 
3.05 - 3,41 
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Position, age and (sub)chron designation of biostratigraphic events according to CK95. Ages have been obtained by linear interpolation 
of the sedimentation rate between dated reversal boundaries (see Appendix A). Ages and (sub)chron designation is given in italics if based 


on cyclostratigraphic correlations to section with magnetostratigraphy. 
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